arXiv: 1501.06403vl [astro-ph.EP] 26 Jan 2015 


Astronomy & Astrophysics manuscript no. TrES-4'letter arxiv © ESO 2015 

January 27, 2015 


Letter to the Editor 


The GAPS Programme with HARPS-N@TNG 
VI: The Curious Case of TrES-4b * 

A. Sozzetti', A.S. Bonomo', K. Biazzo^, L. Mancini^, M. Damasso^ S. Desidera"^, R. Gratton^, A.F. Lanza^, E. 
Poretti^, M. Rainer^, E. Malavolta^, E. Affer^, M. Barbieri^, E.R. Bedin^, C. Boccato^, M. Bonavita^, E. Borsa^, S. 
Ciceri^, R.U. Claudi^, D. Gandolfi^’^, P. Giacobbe', T. Henning^, C. Knapic^, D.W. Eatham^, G. Eodato^^, A. Maggio^, 
J. Maldonado^, E. Marzari'^’^, A.E. Martinez Eiorenzano'^, G. Micela^, E. Molinari^^’^^, C. Mordasini^, V. 
Nascimbeni^, L Pagano^, M. Pedani^^, F. Pepe^"*, G. Piotto'^’^, N. Santos'^’^^, G. Scandariato^, E. Shkolnik^^, and J. 

S on th worth 

(Affiliations can be found after the references) 

Received ??; Accepted ?? 


ABSTRACT 

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. 
A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude K = 51 ± 3 m s“'. The 
derived mass of TrES-4b is found to be Mp = 0.49 ± 0.04 Mj^p, significantly lower than previously reported. Combined with the large radius 
(Rp = 1.84 ;^qpp Rjup) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several 
scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets. 
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1. Introduction 

The class of transiting extrasolar planets (to-date, over 1000 are 
either confirmed or validated) allows for many a study to fur¬ 
ther our understanding of their interiors, atmospheres, and ulti¬ 
mately formation and evolution history (see, e.g. Madhusudhan 
et al. 2014 and Baraffe et al. 2014). The subset of close-in gi¬ 
ant planetary companions (hot Jupiters) with very large radii, 
and corresponding very low mean densities, posed for a time a 
conundrum to theoreticians (the so-called radius anomaly prob¬ 
lem; see e.g. Bodenheimer et al. 2003). There is now a growing 
consensus that the radius of a hot Jupiter can be inflated due to 
several factors, including variable stellar irradiation, the planet’s 
mass and heavy element content, tidal and kinetic heating, and 
Ohmic dissipation (for a review see, e.g., Spiegel et al. 2014, and 
references therein). 

The distribution of planetary radii of transiting hot Jupiters 
in systems with well-determined stellar and planetary param¬ 
eters has been described in the recent past in terms of some 
of the relevant factors (such as equilibrium temperature, stel¬ 
lar metallicity, and orbital semi-major axis) using empirical for¬ 
mulae based on the assumption of independent variables (Beky 
et al. 2011; Enoch et al. 2011) or a multivariate regression ap¬ 
proach (Enoch et al. 2012; Weiss et al. 2013). These latest mod¬ 
els are quite successful in statistically reproducing the observed 


* Based on observations made with the Italian Telescopio Nazionale 
Galileo (TNG) operated on the island of La Palma by the Fundacion 
Galileo Galilei of the INAF at the Spanish Observatorio del Roque 
de los Muchachos of the lAC in the frame of the program Global 
Architecture of Planetary Systems (GAPS), and with the Zeiss 1.23- 
m telescope at the German-Spanish Astronomical Center at Calar Alto, 
Spain 


radius distribution of this class of exoplanets. Still, some of the 
most extreme planets with the largest radii remain challenging 
for current models of planetary formation and bulk structure. 
For example, the extremely low densities of objects such as 
WASP-17b (Anderson et al. 2010), HAT-P-32b (Hartman et al. 
2011), WASP-79b (Smalley et al. 2012), WASP-88b (Delrez et 
al. 2014) or Kepler-12b (Fortney et al. 2011) cannot be repro¬ 
duced by simple models of core-less planets (e.g., Baraffe et al. 
2014), nor can the atmospheric inflation mechanisms mentioned 
above explain the observed radii. 

TrES-4b (Mandushev et al. 2007, M07 thereafter) is another 
highly bloated transiting hot Jupiter. It belongs to the restricted 
lot of some dozen objects with a measured radius larger than 1.7 
J^jup (Sozzetti et al. 2009, S09 thereafter; Chan et al. 2011; Sada 
et al. 2012; Southworth 2012). Measurements of the Rossiter- 
McLaughlin effect (Narita et al. 2010, NIO hereafter) revealed 
close spin-orbit alignment of the TrES-4 system. Atmospheric 
characterization measurements have been obtained by Knutson 
et al. (2009), who detected a temperature inversion in the TrES- 
4b’s broadband infrared emission spectrum with Spitzer/IRAC 
during secondary eclipse, and by Ranjan et al. (2014), who pre¬ 
sented a featureless transmission spectrum of TrES-4b using 
HST/WFC3 during primary transit. Constraints from the sec¬ 
ondary eclipse measurements and expanded radial-velocity (RV) 
datasets (Knutson et al. 2014, K14 thereafter) indicate a proba¬ 
ble circular orbit for TrES-4b. Using a time baseline in excess 
of five years K14 did not detect any significant acceleration in 
the RV data that might point to the presence of a massive outer 
companion in the system. Finally, TrES-4 has a faint common 
proper motion companion at ~ 1.5", discovered by Daemgen et 
al. (2009) and confirmed by Bergfors et al. (2013). 
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Table 1. HARPS-N radial velocities, formal errors, bisector 
spans, and chromospheric activity index of TrES-4. 


BJDtdb 

-2450000 

RV 

(kms^') 

±1 cr 
(kms^') 

Bis. span 
(kms“') 

^ozR'hk 

6362.736988 

-16.069 

0.007 

-0.005 

-5.046 

6484.535732 

-16.139 

0.008 

0.034 

-5.290 

6485.620551 

-16.112 

0.011 

- 0.011 

-5.303 

6506.435512 

-16.143 

0.007 

0.018 

-5.154 

6507.544966 

-16.078 

0.018 

0.050 

-5.181 

6508.534111 

-16.071 

0.010 

-0.045 

-5.259 

6509.530003 

-16.139 

0.009 

- 0.022 

-5.257 

6543.363868 

-16.052 

0.007 

- 0.010 

-5.068 

6583.321741 

-16.036 

0.010 

-0.041 

-5.189 

6586.336543 

-16.039 

0.009 

-0.024 

-5.109 

6696.775058 

-16.078 

0.020 

-0.033 

-5.058 

6699.728820 

-16.049 

0.009 

- 0.012 

-5.117 

6701.762988 

-16.145 

0.009 

0.007 

-5.099 

6786.659718 

-16.135 

0.012 

- 0.010 

-5.160 

6787.693684 

-16.130 

0.010 

0.010 

-5.111 

6857.483512 

-16.109 

0.007 

0.030 

-5.089 

6858.451440 

-16.151 

0.008 

0.027 

-5.169 


In this Letter we present RV measurements of TrES-4 gath¬ 
ered with the HARPS-N spectrograph (Cosentino et al. 2012) 
on the Telescopio Nazionale Galileo (TNG) within the context 
of the programme Global Architecture of Planetary Systems 
(GAPS, Covino et al. 2013; Desidera et al. 2013), along with 
additional photometric light-curves during transit obtained with 
the Zeiss 1.23-m telescope at the German-Spanish Calar Alto 
Observatory (CAHA). A combined analysis allows us to derive a 
much lower mass for TrES-4b than previously reported, making 
it the second-lowest density transiting hot Jupiter known to-date. 

2. Spectroscopic and Photometric Observations 

The TrES-4 system was observed with HARPS-N on 17 indi¬ 
vidual epochs between March 2013 and July 2014. The Th-Ar 
simultaneous calibration was not used to avoid contamination 
of the stellar spectrum by the lamp lines (which would affect 
a proper spectral analysis). In addition, the magnitude of the 
instrumental drift during a night (< 1 m s ') is considerably 
lower than the typical photon-noise RV errors (:^ 9 m s '), thus 
of no impact for faint stars such as TrES-4 (see, e.g., Bonomo 
et al. 2014; Damasso et al. 2015). The reduction of the spectra 
and the RV measurements were obtained using the latest version 
(Nov. 2013) of the HARPS-N instrument data reduction software 
(DRS) pipeline and the G2 mask. We measured the RVs using 
the weighted cross-correlation function (CCE) method (Baranne 
et al. 1996; Pepe et al. 2002). The individual measurements are 
reported in Table [1] together with the values of bisector span and 
chromospheric activity index logR^^. 

The spectra of TrES-4 were coadded to produce a merged 
spectrum with a peak signal-to-noise ratio of ~ 110 pixel ' at 
550 nm. We determined the atmospheric stellar parameters using 
the code MOOG (Sneden 1973; version 2013) and implemented 
both the methods based on equivalent widths and on spectral 
synthesis, as described in Biazzo et al. (2012), D’Orazi et al. 
(2011), and Gandolh et al. (2013). The hnal adopted parameters 
are listed in Table |2] 

We carried out /c-band precision photometric observations 
of two complete transit events of TrES-4 b with the CAHA 1.23- 
m on UT 2013 July 6 and UT 2014 June 30. The telescope was 


Table 2. System parameters of TrES-4. Errors and upper limits 
refer to 1 cr uncertainties. 


Stellar parameters 

Effective temperature TcffiK] 

6295 ± 65 

Metallicity [Fe/H] [dex] 

0.28 ± 0.09 

Microturbulence velocity [ km s“'] 

1.73 ±0.02 

Rotational velocity V sink [kms^'] 

8.5 ± 0.5 

Systemic velocity y [ m s“'] 

-16097.0 ±2.6 

RV jitter [ms^'j 

< 6 

Density p, [g cm“^] 

0 ^ 47+0 042 

U.JJH-/ 0,031 

Mass [Mq] 

1.45 ±0.05 

Radius [Re] 

1.81 ±0.08 

Derived surface gravity log g [cgs] 

4.09 ± 0.03 

Age t [Gyr] 

2.2 ± 0.4 

Transit and orbital parameters 

Orbital period P [days] 

3.55392771 (47) 

Transit epoch TnrBJDxnR - 24500001 

4230.90560 (30) 

e coscri 

u.uuiu_oooi7 

e sin to 

rj+ 0.012 

^- 0.022 

Orbital eccentricity e 

< 0.016 

Argument of periastron at [deg] 

unconstrained 

RV semi-amplitude A" [ms“'] 

51.1 ±3.3 

Transit duration T 14 [h] 

'2 ^co+0.036 

j.ujo_oo30 

Radius ratio Rp/R, 

Inclination i [deg] 

Linear limb-darkening coefficient u 
ajR, 

Impact parameter b 

f) 

u.iu-+jz,_o 00072 
83.071“;^ 
n S 74+0060 

5 J 4 + 0 . 2 T 

0744“ri’o,6 

^•'^^- 0.022 

Planetary parameters 

Mass Mp [Mjup] 

0.494 ± 0.035 

Radius Rp [Rj^p] 

1 

i.ojo_o,o90 

Density Pp [g cm“^] 

u.u:?>'_0.oi3 

Surface gravity log gp [cgs] 

2.45 ± 0.05 

Orbital semi-major axis a [au] 

0.0516 ±0.0005 

Equilibrium temperature [K] 

1795^35 


^ from lMandushev etaP (20071) 

^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat 
redistribution to the night side. 

defocussed and autoguided during all the observations and the 
CCD was windowed to reduce the readout time. The datasets 
were reduced using standard calibration techniques (overscan 
correction, trimming, bias subtraction, flat helding). We then de¬ 
rived differential fluxes relative to an ensemble of local compar¬ 
ison stars (using the methodology described in Southworth et 
al. 2014). The hnal set of photometric time series of TrES-4 is 
available in a machine-readable form in the electronic version of 
Table [3] Uncertainties on individual photometric measurements 
were estimated separately for the two light curves as the stan¬ 
dard deviation of the residuals of the transit htting; these un¬ 
certainties are larger than the formal error bars in both cases. 
Correlated noise was then estimated following Pont et al. (2006) 
and Bonomo et al. (2012) and added in quadrature with the indi¬ 
vidual measurement uncertainties. Einal uncertainties are equal 
to 8.42 X 10“"' and 7.71 x 10“"' (in units of relative hux) for the 
former and latter light curve, respectively. 

3. Revised TrES-4 System Parameters 

New parameters of the TrES-4 system were derived through a 
Bayesian combined analysis of our photometry in U band and 
HARPS-N RV meas u rements by simultaneously htting a tran¬ 
sit model dGimenezl l2006l l2009h and a Keplerian orbit. Eor 
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- 0.04 - 0.02 0.00 0.02 0.04 

Orbital phase 

Fig. 1. Phase-folded observations of two full transits of TrES- 
4b in /c band along with the best-fit model (red solid line). See 
Section 3 for details. 


Table 3. Differential photometry of TrES-4. The full dataset will 
be made available at the CDS. 


4 band (CAHA 1.23-m) 


BJDtdb - 2450000 

Relative flux 

Uncertainty 

6480.353148 . 

1.000482 

0.000842 

6480.355590 . 

1.000129 

0.000842 

6480.357108 . 

0.999889 

0.000842 

6480.358670 . 

0.999212 

0.000842 

6480.360246 . 

1.000970 

0.000842 

6480.361822 . 

1.000361 

0.000842 


this purpose, we use d a Differential Evolution Markov C hain 
Monte Carlo method (iTer Braakl2006tlEastman e t al.l2013T) with 
a Gaussian likelihood function (see, e.g., lGregorvll2005h . Our 
global model has eleven free parameters: the transit epoch Tq\ 
the orbital period P; the systemic radial velocity y; the radial- 
velocity semi-amplitude K\ e cos u and e sin w, where e is the ec¬ 
centricity and CO the argument of periastron; an error term added 
in quadrature to the formal uncertainties to account for possi¬ 
ble jitter in the RV measurements regardless of its origin, such 
as instrumental effects, stellar activity, etc.; the transit duration 
from first to fourth contact T 14 ; the ratio of the planet to stellar 
radii Pp/P*; the inclination i between the orbital plane and the 
plane of the sky; and the coefficient m of a linear limb-darkening 
law. We first tried to use a quadratic limb-darkening law but the 
two coefficients, especially the quadratic one, were highly un¬ 
constrained. This means that the precision of our transit light 
curves does not allow to fit both coefficients. 


Gaussian priors were imposed on Pq and P, after improving 
the transit ephemeris by combining the transit epochs available 
in the literature (M07; Chan et al. 2011) with the two epochs 
we derived from our 4 photometry by analyzing each individ¬ 
ual transit with a circular transit model and a DE-MCMC tech¬ 
nique. Gaussian priors were als o set on the center tim es of the 
secondary eclipses observed by iKnutson et alJ (l2009h with the 
Spitzer space tel escope because these p rovide strong constraints 
on e cos to (e.g.. [Jordan & Bakosll2008h . Non-informative priors 



Orbital phase 

Fig. 2. Top panel: Phase-folded RV measurements of TrES- 
4 obtained with HARPS-N (blue circles) and, superimposed, 
the best-fit Keplerian orbit model (black solid line). The 
Keck/HIRES RVs from M07 (green diamonds) and K14 (red 
squares) and the two best-fit orbital solutions obtained in those 
papers are also shown. Bottom panel: Residuals from the best-fit 
model to the HARPS-N radial velocities. 


were used for the other orbital and transit parameters while a 
modified Jeffrey’s prior was adopted for the RV jitter term. 

The DE-MCMC analysis was stoppe d afte r reach ing conver¬ 
gence and good mixing of the chains dEordl l2006h . The final 
best-fit transit model and RV curve are overplotted on the phase- 
folded data in Eigures [T] and |2l respectively. We do not deter¬ 
mine a significant RV jitter, listed as an upper limit in Table |2] 
(only internal errors are reported in Table [1] and Eigure|2]i. The 
density of the host star from the transit fitting, and the effec¬ 
tive temperature and stellar metallicity as derived in Sect. 2 
were later compared with the theo retical Yonsei-Yale evolution¬ 
ary tracks (iDemaraue et al.L [200411 to determine the stellar mass, 
radius, surface gravity, age, and t h eir as sociated uncertainties 
(ISozzetti et al.L |2007[ iToires et’all l2012h . These are listed in 
Table |2l and agree within Icr with th e literature values (cf., e.g., 
iTorre^t"^ l2()08[ IChan et ^1201 ih . The related planetary pa¬ 
rameters are Mp - 0.494 + 0.035 Mjup, Rp - 1.838j;QQgy Rjup, 
and pp - 0.099;(Q y [3 g cm“^. 


4. Discussion and Conciusions 

Our new, combined spectroscopic and photometric analysis of 
the TrES-4 system allows us to determine stellar properties in 
good agreement (within the errors) with those measured by M07 
and S09. The planetary radius agrees also well with the most 
recent determinations by S09 and Chan et al. (2011), although 
we formally derive its largest value to-date. However, one strik¬ 
ing element emerges from our study. The best-fit Keplerian orbit 
for TrES-4 based on HARPS-N RV measurements (K = 51 + 3 
m s ') has an amplitude almost a factor 2 smaller than the one 
{K - 91 ±1 rc\ s“') reported by M07. As a consequence, the re¬ 
vised mass of the planet is ~ 1.7 times lower. Such a discrepancy 
clearly deserves a thorough investigation, and we describe here 
the steps we have taken in this direction. 
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In the most recent update of the TrES-4 system parameters, 
K14 report K - 84 + 10 m s“* and Mp = 0.843^°°^^ Mjup. These 
numbers are compatible within the error-bars wiih the initial es¬ 
timates of M07 and S09. The K14 orbital solution is based on 
the combination of three datasets, from M07, NIO, and obtained 
by the authors themselves. We show in Figure |2] a phase plot of 
the published Keck velocities and our HARPS-N dataset, super¬ 
posed on the three orbital solutions derived by M07, K14, and in 
this work (the NIO RV set obtained with Subaru/HDS is of sig¬ 
nificantly lower internal precision and is not shown). The RVs 
published by M07 are clearly incompatible with the K-value de¬ 
rived based on HARPS-N RV data. The data obtained by K14 
did not sample the critical orbital phases, and they appear con¬ 
sistent with both solutions. The larger K-value in the Keck data 
is thus driven by the observations of the discovery paper. 


2.0 
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Fig. 3. Mass-radius diagram of the known transiting planets with 
Rp > 0.4 Rjup and Mp >0.1 Mjup. Only systems with masses de¬ 
termined to better than 30% precision are included. Green dia¬ 
monds indicate the Solar System giant planets Jupiter and Saturn 
(from right to left). The three dotted lines display isodensity 
curves of 0.1, 0.5, and 1.5 gcm“^. The position of TrES-4b as 
determined in this work is shown with a hlled red square, to be 
compared with that (empty red square) derived by S09. 


There are several scenarios that can be proposed to explain 
the observed discrepancy in the RV amplitudes. One possible 
culprit might be the faint companion at ~ 1.5" (almost due north 
of TrES-4). Cunha et al. (2013) have analyzed in detail the im¬ 
pact on precision RVs of faint stellar companions from spectra 
gathered with fiber-fed spectrographs. The companion of TrES- 
4 is of late-K or early-M spectral type and 4.5 mag fainter 
at /-band. From Table 8 of Cunha et al. (2013) one then infers 
that contamination levels between 1 and 10 m s * could apply 
in case the companion were to fall within the hber of HARPS- 
N. A systematic effect of similar magnitude might be induced 
on the Keck RVs if the companion had been on the HIRES slit 
during the period of the M07 observations. Either way, this sce¬ 
nario does not seem to provide a convincing explanation for the 
observed difference in the K-value, as the higher RV dispersion 
inferred does not have the required magnitude and, most impor¬ 
tantly, such effect would have had to occur in such a way as to 
exactly double (or halve) the signal amplitude. The hypothesis 
of large starspots on the stellar photosphere (causing an apparent 
RV shift on a timescale of the stellar rotation period) is unlikely 
for a late F-star such as TrES-4 (see Knutson et al. 2009). Large, 
intrinsic stellar jitter also does not appear to be supported by the 
observational evidence. No emission is seen in the Can H&K 
lines related to magnetic activity in the HARPS-N spectra, from 



Fig. 4. The dependence on equilibrium temperature of the ob¬ 
served radii of giant planets with Mp > 0.1 Mjup and Rp > 
0.8 Rjup (data from http://exoplanet.eu/). Purple squares, grey 
circles, and green triangles indicate objects with pp < 0.25 g 
cm~^, 0.25 < Pp < 1.5 g cm“^, and pp > 1.5 g cm“^, respec¬ 
tively. The location of TrES-4b is shown with a red square. 


which we derive (logRj^^) = -5.15 + 0.08, essentially indistin¬ 
guishable from the value reported by S09. The empirical relation 
by Wright (2005) predicts a typical stellar jitter of ~ 4 m s“' for 
such a star as TrES-4. We note however that the very low value 
of the chromospheric emission is in line with the correlation 
found by Hartman (2010) and could be explained as the effect 
of absorption in the Ca ii H&K line cores by material evaporated 
from the low-gravity planet (Lanza 2014; Figueira et al. 2014). 
Based on the absence of bump progression in the mean line pro¬ 
files of HARPS-N and archival Keck data (determined with the 
Donati et al. 1997 technique), we also ruled out the possibility 
that the amplitude of the RV curve of TrES-4 would be modified 
by non-radial stellar pulsations typical of y Dor variables (Kaye 
et al. 1999), which have been detected in a few cases in stars 
with similar stellar parameters to those of TrES-4 (Uytterhoeven 
et al. 2014). A homogeneous, comprehensive re-analysis of all 
available Keck data on the system might help to resolve the co¬ 
nundrum, particularly to understand if unrecognized systematics 
in the few Keck RVs in the discovery paper could be to blame. 

The much smaller mass of TrES-4b as determined by 
HARPS-N implies a significantly lower density than previ¬ 
ously thought for this hot Jupiter. The new location of TrES- 
4b in the mass-radius diagram of known transiting giant planets 
(Figure|3]) makes it the second-lowest density object, WASP-17b 
(Southworth et al. 2012; Bento et al. 2013) being the record- 
holder at present. With a mass closer to Saturn’s, TrES-4b’s 
predicted radius is signihcantly underestimated by all empirical 
relations recently proposed in the literature (Beky et al. 2011; 
Enoch et al. 2011, 2012; Weiss et al. 2013), with radius differ¬ 
ences ranging between 0.74 Rjup and 0.45 Rjup. In a Rp - Tgq di¬ 
agram (see Figure|4]for details) TrES-4b nicely hts in the upper 
envelope of lowest-density objects, which exhibit a strong pos¬ 
itive correlation between the two parameters (Spearman’s rank 
correlation coefficient = 0.82 + 0.03). We note how the trend 
of increasing Rp with Teq becomes significantly milder if we con¬ 
sider the sample of the densest giants (r^ - 0.54 + 0.06), and the 
relationship becomes completely flat if a cut-off around 1.0 Rjup 
is adopted (rather than the one used in Figure|4]). 

We conhrm a very low eccentricity (e < 0.016 at the Icr 
level) for TrES-4b’s orbit, improving upon the recent determi¬ 
nation by K14. An estimate of the typical tidal timescales based 
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on the model by Leconte et al. (2010) adapted as to allow con¬ 
stant modified tidal quality factors for the star (2' = 10^) and 
the planet {Q'^ - 10^) gives a circularization timescale of 40 Myr 
supporting the e 0 hypothesis. On the other hand, the timescale 
for the evolution of the obliquity obtained with the same tidal 
model is ~ 24 Gyr, while that for the orbital decay is ~ 6 Gyr. 
This suggests that the alignment of the system is primordial and 
that no remarkable tidal evolution of the orbit has occurred dur¬ 
ing the main-sequence lifetime of the system. With the presently 
derived upper limit for the eccentricity, the maximum power dis¬ 
sipated by equilibrium tides inside the planet is ~ 4.5 x lO'^ W, 
insufficient to explain its large radius anomaly. Given its pecu¬ 
liarity, further photometric and spectroscopic monitoring of the 
TrES-4 planetary system is clearly encouraged. 
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